Chiral properties of strong interactions in a magnetic background 
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We investigate the chiral properties of QCD in presence of a magnetic background field and in 
the low temperature regime, by lattice numerical simulations of Nf = 2 QCD. We adopt a standard 
staggered discretization, with a pion mass around 200 MeV, and explore a range of magnetic fields 
(180 MeV) 2 < \e\B < (700 MeV) 2 , in which we study magnetic catalysis, i.e. the increase of chiral 
symmetry breaking induced by the background field. We determine the dependence of the chiral 
condensate on the external field, compare our results with existing model predictions and show that 
a substantial contribution to magnetic catalysis comes from the modified distribution of non-Abelian 
gauge fields, induced by the magnetic field via dynamical quark loop effects. 

PACS numbers: 12.38. Aw, 11.15.Ha,12.38.Gc 



I. INTRODUCTION 

The study of strong interactions in presence of a strong 
background magnetic field has attracted increasing at- 
tention in the recent past. On one side the issue is 
of great phenomenological relevance: magnetic fields 
of the order of 10 16 Tesla, i.e. y/\^\B ~ 1.5 GeV) 
may have been produced at the cosmological electrowcak 
phase transition [l| and they may have influenced sub- 
sequent strong interaction dynamics, including the con- 
finement /deconfmement transition. Slightly lower fields 
are expected to be produced in non-central heavy ion 
collisions (up to 10 14 Tesla at RHIC and up to - 10 15 
Tesla at LHC @, 0]), where they may give rise to new 
phenomenology, the so-called "chiral magnetic effect", 
capable of revealing the presence of deconfined matter 
and of non-trivial topological vacuum fluctuations ■ 
Finally, magnetic fields of the order of 10 10 Tesla are ex- 
pected to bepresent in a class of neutron stars known as 
magnetars Q (for a recent review see Ref. 0). 

On the other side, a background magnetic field 
(electro-magnetic or chromo-magnetic) may serve as yet 
another parameter to probe the structure of the QCD 
vacuum and of the QCD phase diagram, on the same 
footing with other interesting external conditions such as 
a baryon chemical potential. Many studies [§-[34j have 
investigated the chiral properties of the theory and what 
is generally known as magnetic catalysis, consisting in an 
enhancement of chiral symmetry breaking and in sponta- 
neous mass generation induced by the magnetic field, a 
phenomenon predicted from different low energy models 
and approximations of QCD and related to the dimen- 
sional reduction taking place in the dynamics of particles 
moving in a strong external magnetic field [13, [H, HH ■ 
More recently, the issue of the influence of a magnetic 
field on the deconfmement transition has been investi- 
gated by means of both lattice QCD simulations and low 
energy models of strong interactions [H], l33 - l47| : there 



is converging evidence that a magnetic field leads to an 
increase of both the strength and the temperature of 
the transition; in the case of a chromo-magnetic field, 
instead, numerical simulations show a decrease of the 
transition temperature [48|,|4!|. Finally, conjectures have 
been proposed according to which a strong enough mag- 
netic field may induce the appearance of new supercon- 
ductive phases [5fj| - l52| . 

In the present study we address the issue of magnetic 
catalysis, presenting the first study of such phenomenon 
by lattice QCD simulations which include the contribu- 
tion of dynamical quarks; previous numerical studies in- 
deed have only considered the effect of the magnetic field 
on quenched configurations [H, H3|- I n particular, we 
have considered QCD at zero or low temperature, with 
two dynamical flavors carrying different electric charges, 
corresponding respectively to the u and d quark charges, 
and coupled to a background constant and uniform mag- 
netic field. We have adopted a standard rooted staggered 
fermion discretization, with a (Goldstone) pion mass of 
about 200 MeV. 

One of the purposes of our investigation is to obtain in- 
formation about the dependence of the chiral condensate 
on the magnetic field, and compare it with various exist- 
ing model and low energy predictions. A second purpose 
that we have is to understand which part of magnetic 
catalysis is a purely tree level effect, due to the fact that 
quarks propagate in a modified background obtained by 
adding the U(l) field to the non-Abelian gauge configura- 
tions, and which part is due to a modification of the non- 
Abelian fields themselves, induced by the loop effects of 
dynamical quarks coupled to the magnetic background. 
Both effects can in principle modify the spectrum of the 
Dirac operator, leading to a larger density of eigenvalues 
around zero, hence to an increase of the chiral condensate 
via the Banks - Casher relation [55| . 

To that aim, one could compare with existing in- 
vestigations of magnetic catalysis, based on SU(2) and 
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SU(3) gauge configurations sampled by pure gauge sim- 
ulations [H, [H| : however that would not be completely 
satisfactory and would also be difficult because of differ- 
ent scale settings and renormalization effects. What we 
will do instead is to try separating the two different ef- 
fects, by inserting alternatively the magnetic field only 
in the computation of the quark propagator (sampling in 
this case configurations without the presence of the mag- 
netic field), or only in the sampling measure, i.e. in the 
fermion determinant, without affecting the quark propa- 
gator computation. We will call "valence" catalysis the 
first contribution and "dynamical" catalysis the second: 
both of them will be compared with the full increase of 
the chiral condensate, obtained when the magnetic field 
is inserted directly both in the fermion determinant and 
in the computation of the quark propagator. As we will 
show, the purely dynamical contribution corresponds to 
a considerable part of the total increase in the quark con- 
densate. 

The paper is organized as follows. In Sec. [TT] we give 
some details about our lattice discretization of QCD in 
presence of a magnetic field and about our numerical 
setup. In Section IIIII we present our numerical results 
and finally, in Section [IV! we give our conclusions. 



II. NUMERICAL SETUP 

The discretization of Nf = 2 QCD in presence of a 
magnetic background field adopted in the present work 
is similar to that reported in Ref. [42| . In particular, 
partition function of the (rooted) staggered fermion dis- 
cretized version of the theory in presence of a non-trivial 
electro-magnetic (e.m.) background field and with differ- 
ent electric charges for the two flavors, q u = 2|e|/3 and 
q c i = — |e|/3 (|e| being the elementary charge), is written 
as: 

Z(T, B) = J VUe~ Sa det M 1 [B, q u ] det M* [B, q d ] (1) 



^ 4 

Mij[B, q] = amSij + — (u(B, q)i <u Ui tV 8i t j-b 

"'I I !■<!>, ; „r ,,A ,.,) • (2) 

T>U is the functional integration over the non-Abelian 
gauge link variables J7 n , M , Sg is the discretized pure 
gauge action (we consider a standard Wilson action); 
u(B, q)i M are instead the Abelian gauge links correspond- 
ing to the background e.m. field. The subscripts i and 
j refer to lattice sites, v is a unit vector on the lattice 
and Tji )V are the staggered phases. Periodic (antiperiodic) 
boundary conditions (b.c.) must be taken, in the finite 
temperature theory, for gauge (fermion) fields along the 
Euclidean time direction, while spatial periodic b.c. are 
chosen for all fields. 



We shall consider a constant and uniform magnetic 
field B = Bz. The presence of periodic b.c. in the x and 
y directions imposes a constraint on the admissible values 
of B, which get quantized, as illustrated in the following 
subsection. Symmetry under charge conjugation imposes 
that Z(T, B) as well as other charge even observables, 
including the chiral condensate, be even functions of B. 



A. Magnetic Field on a Torus 

In presence of periodic b.c, the magnetic field in the 
z direction goes through the surface of a torus in the 
x — y directions, whose total extent is l x l y - The circu- 
lation of along any closed path, lying in the x — y 
plane and enclosing an arbitrary region of area A (see 
e.g. Fig. [T]), is proportional, by Stokes' theorem, to the 
flux of B through the enclosed surface 



A^dx^ = AB 



(3) 



On the other hand, since we are on a torus, it is ambigu- 
ous to state which is the enclosed surface: the comple- 
mentary region of area l x l y — A can be chosen as well, 
therefore one can equally state 

<j> A^dx^ = (A — lJ y )B . (4) 

At the level of the gauge field the ambiguity is resolved 




FIG. 1: Surface with periodic b.c. orthogonal to the direction 
of the magnetic field. The phase factor taken by a particle 
moving around the plotted countour must be defined unam- 
biguously and this leads to magnetic flux quantization. 

by admitting discontinuities in A^ somewhere on the 
torus, or alternatively by covering the torus with various 
patches where different gauge choices are taken. In any 
case, one has to guarantee that the ambiguity is not visi- 
ble by charged particles moving on the torus, and this is 
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true only if the phase factor taken by the charged particle 
moving along the closed path is defined unambiguously 



i.e. if 



exp (iqBA) = exp (iqB(A — l x ly)) 



qB = 2irb/l x l y 



(5) 



(0) 



where b is an integer. Notice that this line of argument is 
exactly the same that applies on a sphere and that leads 
to Dirac quantization of the magnetic monopole charge. 
The quantization rule depends on the electric charge of 
the particles feeling the presence of the magnetic field, in 
particular it is set by the smallest charge unit, which in 
our case is brought by the d quark, q d = — |e|/3, hence 



\e\B = 6nb/l x l y = 6nba 2 / L x L y 



(7) 



where L x and L y are the system sizes in lattice units and 
a is the lattice spacing. 

Further details about the definition of a magnetic field 
on a torus can be found in Ref. (59j , where it is shown that 
translational invariance on the torus is explicitly broken 
by the presence of the magnetic field, due to the non- 
trivial phases taken by particles winding around one of 
the directions of the torus (Wilson lines) : only a discrete 
invariance is left, by shifts which are integer multiples of 



= l x /b 



ly/b 



(8) 



respectively in the x and y directions. Such invariance 
is reduced further on a lattice, since only shifts, if any, 
which are multiples of both a x (a y ) and the lattice spac- 
ing a leave the system invariant: that may lead to addi- 
tional discretization effects. 



B. Discretization details 

We have taken the following choice for the continuum 
e.m. gauge field: 

A y = Bx ; A^ = for ji = x, z,t . (9) 

The corresponding U(l) links on the lattice are: 

u{B,q) n , y = e la ' 2qBn * ■ u(B,q) n ^ = 1 for fj. = x,z,t {10) 

In order to guarantee the smoothness of the background 
field across the boundary and the gauge invariance of the 
fermion action, the U(l) gauge fields must be modified 
at the boundary of the x direction: 



u(B,q) n 



. — ia qL x Bn y 

x \n x — L x ~ " 



(ii) 



and the magnetic field must be quantized as specified 
in Eq. ||7J). That corresponds to taking the appropriate 
gauge invariant b.c. for fermion fields on the torus 
(with the possible additional free phases 9 X and 9 y 
set to zero). 



59] 
59] 



We have considered a symmetric lattice, L x = L y = 
L z = N t = 16, a bare quark mass am — 0.01335 and 
an inverse gauge coupling /3 = 5.30. According to scale 
estimates reported in Ref. [42] ]. that corresponds to a 
lattice spacing a ~ 0.3 fm, a (Goldstone) pion mass m, ~ 
200 MeV and a temperature T = (i^a)" 1 ~ 40 MeV, 
hence low enough that the system can be considered to 
be effectively at zero temperature. 

We have explored different values of \e\B which, ac- 
cording to Eq. ([7]), can be changed only in units of 
6ira~ 2 / L x L y ~ (180 MeV) 2 . Notice however that the 
presence of an ultra-violet (UV) cutoff imposes also an 
upper limit on the possible values of B that can be ex- 
plored on the lattice. To appreciate that, let us consider 
again the phase factor picked up by a particle moving 
around a closed path in the x — y plane, and which con- 
tains all the relevant information about the effect of the 
magnetic field on particle dynamics: there is a minimal 
such path on the lattice, corresponding to a plaquette, 
around which the particle takes the phase factor 



exp(iqa 2 B) = exp 



ilub 

L X Ly 



(12) 



The phase factor above, and therefore all other phase fac- 
tors associated to any closed lattice path, cannot distin- 
guish magnetic fields such that qa 2 B differs by multiples 
of 2ir. One can therefore define a sort of "first Brillouin 
zone" for the magnetic field, 



i.e. 



7T 7T 

— <qB < — 



L x L y ^ < L x L y 



(13) 



(14) 



with all physical quantities being periodic in qB (b) with 
a period 1*k j a 2 (L x L y ); symmetry under b — > —b further 
reduces the range of interesting values of b. Even before 
reaching the limits reported in Eq. (|13[) . one expects the 
periodicity to induce saturation effects, which may dis- 
tort the true physical dependence of observables, like the 
chiral condensate, on B. One should always worry about 
the possible presence of such saturation effects, when try- 
ing to extract information relevant to continuum physics. 



C. Observables and simulation details 

The quantity which is the subject of our investigation 
is the chiral condensate. In presence of a non-zero B wc 
can define two different condensates 



j/d 



(B) 



d\ogZ t 

dm u/d 



\m u/d =m 



(15) 



= yp[/7'[m,C/,i?]Tr(M- 1 [^S,'? u /<i] 
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where the functional integral measure is (see Eq. ([T])): 
V [to, U, B] cx det MJ[m, B, q u ]det M *[m,B, q d ]e- s ° (16) 

A quantity which is useful to discuss magnetic catalysis 
is the relative increment of the quark condensate, which 
we define as: 

E u/d (B) S u/d (B) - S(0) 
r u /d{B) = — ^ 1 = . (17) 



E(0) 



E(0) 



The advantage of r(B) is that it is a dimensionless quan- 
tity and that most renormalizations appearing in the def- 
inition of E cancels out in Eq. (|T7|) . Indeed, assuming 
that renormalizations have a negligible dependence on 
B, as should be the case as long as B stays away from 
the scale of the UV cutoff, the mass dependent additive 
renormalization of E will cancel out in the numerator of 
Eq. (1171) . A residual additive renormalization remains in 
the denominator, leading to an incorrect overall normal- 
ization of r(B): in the following we shall try to estimate 
the magnitude of such systematic error. 

We shall define flavor averaged quantities as well: 



E(B) = 



Z u (B)+Z d (B) 



and 



r(B) 



gCg) 

E(0) 



r u {B)+r d (B) 



(18) 



(19) 



Both E U (_B) and T, d (B) are, by charge conjugation 
symmetry, even functions of B. Moreover, according to 
what discussed in Sec. Ill Bt they are expected to be peri- 
odic in B, with a period 2ira~ 2 /qd (or alternatively with 
a period L x L y in terms of the quantum number b defined 
in Eq. ©). One could expect the u quark condensate to 
have a periodicity shorter by a factor 2, since \q u \ — 2\q d \, 
however this is not exactly true because of the measure 
P[m, U, B] appearing in Eq. (ITS"]) , whose periodicity is set 
by the quark with the lower charge. 

In the limit to — > the chiral condensate is an order 
parameter for chiral symmetry breaking and is related 
by the Banks - Casher relation [55[ to the density p(X) of 
eigenvalues of the Dirac operator, D — M — to Id, around 
A = 0: E = 7rp(0). On the contrary, for m ^ 0, the chi- 
ral condensates defined in Eq. (TTSl) are not related to the 
densities of zero eigenvalues of the respective Dirac oper- 
ators, p u /d(0), which instead can be obtained by taking 
the limit to — ¥ only for the trace term appearing in 
Eq. (JTSJ), i.e. 



Pu/d(0) = - lim fvUV[m,U}Ti(M~ 1 [m',q u/d })(20) 

where the dependence on B has been left implicit. Wc 
shall consider also such quantities and the corresponding 
relative increments 



As discussed in the introduction, we are also inter- 
ested in studying contributions to magnetic catalysis 
coming separately either from the change in the observ- 
able Tr (M _1 [m, B, q u / d }) ("valence" contribution), or 
from that in the measure V[m, U, B] ("dynamical" con- 
tribution). For that reason we define also: 



T>UV[m, U, 0]Tr (M"V, B, q u/d }) (22) 



and 



E*£(fl) = J VUV[m,U,B]Ti(M- 1 [m,Q,q u /d}) .(23) 

In the first case we look at the spectrum of the fermion 
matrix which includes the magnetic field explicitly, but is 
defined on non-Abelian configurations sampled at B = 0. 
In the second case we look at the spectrum of the fermion 
matrix without an explicit magnetic field, but defined 
on gauge configurations sampled in presence of the mag- 
netic field. From E™' d (_B) and E^^(_B) we can define 

the corresponding quantities, Yl val / dyn , ™, r val / d v n ^ 

analogously to what done in Eqs. ([XT]). (fTgj) and p^|) . 

On general grounds we may expect that, in the 
limit of small fields, B acts as a perturbation for 
both the measure term V[m,U, B] and the observable 
Tr [M- X [m, B, q u/d \) in Eq. {T2J}. Given that both func- 
tions are even in B and assuming they are also analytic 
(this may not be true in some limits, see discussion be- 
low), so that the first non-trivial term in B is quadratic, 
one can write, configuration by configuration: 



V[m,U,B]=V[m,U,0] + C B 2 + 0(B*) (24) 



and 



Tr (Ar 1 ^}) = Tr (M -1 [0]) + C B 2 + 0{B A ) (25) 

where it is assumed implicitly that the two constants C 
and C depend on the quark mass and on the chosen 
configuration. Putting together the two expansions one 
obtains 



E(0) 



1 _ r val 
~ 1 — r u/d 



(B) 



dyn 
u/d 



(B)+0{B 4 ). (26) 



r u/d = Pu/d(0,B)/p(0,B = 0)-l. 



(21) 



Therefore, at least in the limit of small fields, the sep- 
aration of magnetic catalysis in a valence part and in 
a dynamical part is a well defined concept. As we will 
show in the following, that continues to be true, within a 
good approximation, for a large range of fields explored 
in the present study. Notice that an approximate addi- 
tivity of r val and r dyn , like in Eq. (|26|) . would be true 
also for different small field dependences, e.g. linear, in 
Eq. (l24l) and Eq. (|25l) ; hence the assumption above is 
stronger and also implies that magnetic catalysis should 
be a quadratic effect in B, at least for small fields and if 
the partition function is analytic in B = 0. 
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b 


r u {b) 


r d Q>) 


< al {b) 


rT'(b) 


<TJ b ) 


1 


0.0005(18) 


-0.001(2) 


0.0017(20) 


0.0008(20) 


-0.001(2) 


2 


0.0077(19) 


0.0022(20) 


0.0070(19) 


0.0027(18) 


0.0003(21) 


3 


0.0202(16) 


0.0077(18) 


0.0151(20) 


0.0037(19) 


0.0046(21) 


4 


0.0356(22) 


0.0162(23) 


0.0266(19) 


0.0052(19) 


0.0097(23) 


5 


0.0567(18) 


0.0274(20) 


0.0407(19) 


0.0121(19) 


0.0162(26) 


6 


0.0760(19) 


0.0358(20) 


0.0579(20) 


0.0165(18) 


0.0182(23) 


7 


0.0996(16) 


0.0481(16) 


0.0759(20) 


0.0217(19) 


0.0273(18) 


8 


0.1246(17) 


0.0613(18) 


0.0949(19) 


0.0281(19) 


0.0361(20) 


9 


0.1474(16) 


0.0717(18) 


0.1144(19) 


0.0352(18) 


0.0413(18) 


10 


0.1736(17) 


0.0864(17) 


0.1340(19) 


0.0412(18) 


0.0470(19) 


11 


0.2005(18) 


0.1021(18) 


0.1554(19) 


0.0503(19) 


0.0594(23) 


12 


0.2258(16) 


0.1173(16) 


0.1765(19) 


0.0584(19) 


0.0655(19) 


13 


0.2501(17) 


0.1312(17) 


0.1983(20) 


0.0676(18) 


0.0733(22) 


14 


0.2737(18) 


0.1450(17) 


0.2192(20) 


0.0762(20) 


0.0802(25) 


16 


0.3227(19) 


0.1769(18) 


0.2568(20) 


0.0957(19) 


0.0971(21) 


24 


0.4636(23) 


0.2830(25) 


0.3809(21) 


0.1777(19) 


0.1399(34) 


32 


0.5462(22) 


0.3727(22) 


0.4472(24) 


0.2594(21) 


0.1722(28) 


48 


0.6485(22) 


0.5053(22) 


0.5308(23) 


0.3816(21) 


0.2027(28) 


64 


0.6855(23) 


0.5790(23) 


0.5652(24) 


0.4460(21) 


0.2199(30) 


80 


0.6545(22) 


0.6198(23) 


0.5317(23) 


0.4924(22) 


0.2159(28) 


96 


0.5726(21) 


0.6504(22) 


0.4480(22) 


0.5297(22) 


0.2128(26) 


112 


0.3868(19) 


0.6589(22) 


0.2603(21) 


0.5549(23) 


0.1876(22) 


128 


0.1333(18) 


0.6376(22) 


0.0000(19) 


0.5642(22) 


0.1358(25) 


144 


0.3828(21) 


0.6567(22) 


0.2583(21) 


0.5558(22) 


0.1848(23) 



TABLE I: Relative increment of the u and d quark conden- 
sates for various magnetic field values. We report full data, 
as well as valence and dynamical contributions separately. 



We have made use of a Rational Hybrid Monte- Carlo 
algorithm to simulate rooted staggered fermions. Typi- 
cal statistics are of the order of 3k thermalized molecular 
dynamics trajectories for each value of the magnetic field. 
The trace of the inverse of the fermion matrix, appearing 
in Eq. (|15p. has been computed, for each quark flavor, by 
means of a noisy estimator, extracting 10 random vec- 
tores for each configuration and for each value of the 
parameters. Numerical simulations have been performed 
on the apeNEXT facilities in Rome. 



B = and different values of the quark mass, keeping 
the lattice size unchanged, in order to determine the de- 
pendence of £(£? = 0,m) on m; results are reported in 
Fig. O The expected leading order dependence on m in 
the chirally broken phase is the following (see e.g. the 
discussion in Ref. [60j): 

E(S = 0, m) = S(0, 0) + cysVm + c x m + 0(m 3 ) (27) 

where the non-analytic square root term is expected 
from the presence of Goldstone mode fluctuations 
IBH, while the leading order, quadratically divergent con- 
tribution to the additive renormalization affects the lin- 
ear term in m. A fit according to Eq. (|27| gives T,(B = 

0. 0) = 0.565(6), ci/2 = 0.61(9) and a = 2.9(4), with 
X 2 /d.o.f. = 5.7/7, from which we infer that the linear 
term in m accounts for about 6% of the total signal mea- 
sured at the quark mass explored in our investigation, 

1. e. am = 0.01335. 

We conclude that our determinations of r, r val and 
r dyn are digtorted by a common and B independent over- 
all normalization factor, which leads to a systematic ef- 
fect of the order of 10% and does not affect issues such 
as the separation of magnetic catalysis into a valence and 
dynamical contribution, as discussed later in this Section. 



0.75 



0.7 



0.65 



0.6 



0.55 



0.005 0.01 



0.015 0.02 
a m 



0.025 0.03 



III. NUMERICAL RESULTS 

We report in Table U the relative increments of the u 
and d quark condensates respectively, including also mea- 
surements of the valence and of the dynamical contribu- 
tion. Notice that r dyn (B) is exactly the same, by defini- 
tion, for u and d quarks, since in this case the magnetic 
field affects only the fermion determinants. One can also 
explicitly verify from the table that r v u al (B/2) = r v d al {B) 
within errors: this is expected since \q u \ = 2\q ( i\ and since 
in this case the sampling measure is independent of B. 

As stressed in Section Hi CI all definitions of r are af- 
fected by a systematic overall normalization factor, due 
to an additive, mass dependent renormalization present 
in the condensate computed at zero magnetic field, E(0) 
(see Eq. (fTT)) ). In order to estimate the magnitude of 
such systematic effect, we have performed simulations at 



FIG. 2: Dependence of the chiral condensate on the bare 
quark mass at B = 0, together with a best fit curve according 
to Eq. ([27) . 



A. Periodicity and saturation effects 

In Fig. [3] we report results for the normalized con- 
densates T, u / d {B)/T 1 {B) (i.e. 1 + r u / d {B)) over the 
whole range of possible independent values of £?, i.e. 
for b/(L x L y ) ranging from to 1 (see discussion at 
the end of Sec. Ill Bp . Notice that data reported for 
b/(L x L y ) > 0.625 are not the result of direct simula- 
tions, but have been obtained by enforcing the expected 
invariances under b — > b+L x L y (i.e. the above mentioned 
periodicity) and under b — > —b, which put together mean 
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invariance under b/(L x L y ) — > 1 — b/(L x L y ). However, 
such invariance has been verified explicitly for a couple 
of points, b/(L x L y ) = 0.4375 and b/(L x L y ) = 0.5625, for 
which independent simulations have been performed. 

Saturation effects, which are present for large values of 
B, are clearly visible from Fig. [31 where we have also 
reported for comparison the results of two fits to the 
small B region. In particular, we infer from the figure 
that one should keep b/(L x L y ) well below 0.1 in order 
that such effects stay negligible: that means \e\B below 
1/a 2 ~ (700 MeV) 2 in our case. In the following only 
data obtained for b < 16 will be considered as reasonably 
free of saturation effects. 

Notice that the u quark condensate shows an approxi- 
mate periodicity in B which is halved with respect to the 
d quark. That comes from the fact that \q u \ = 2|g<j| and 
is only approximate since instead the measure term has 
the usual periodicity. For instance r u (b) has a minimum 
but is not exactly zero at b/(L x L y ) — 1/2, where the ef- 
fective magnetic field felt by u quarks is zero: there is a 
residual catalysis induced by dynamical d quarks, which 
instead at b/{L x L y ) = 1/2 feel the maximum possible 
magnetic field; for the same reason r<j(&) does not reach 
its maximum at b/(L x L y ) = 1/2. Such effects, which 
are absent for the purely valence contribution as can be 
checked from Table U are a first example of the dynami- 
cal contribution to magnetic catalysis, which we discuss 
in detail in the next subsection. 




FIG. 3: Normalized u and d quark condensates as a function 
of the magnetic field for the whole range of independent possi- 
ble values of B. Data for b/(L x L y ) > 0.6 have been obtained 
by enforcing the expected symmetry of the chiral condensate 



under b/(L x L y 



1 — b/(L x L y ), while such symmetry has 



been verified for a couple of points, b/(L x L y ) = 0.4375 and 



b/(L x L y 



0.5625. We also report two curves correspond- 



ing to best fits in the small field region, to better show the 
presence of saturation effects. 



B. Dynamical and Valence contributions to 
magnetic catalysis 

In Fig.Hwe report the functions r(B), r val (B), r d y ,l (B) 
(see Eqs. and {23$) as well as the sum r val {B) + 
r dyn (B), in order to appreciate the amount of magnetic 
catalysis caused by the modified distribution of the non- 
Abelian gauge fields, induced by the coupling of dynam- 
ical quarks to the magnetic field (dynamical contribu- 
tion). We have limited our analysis to b < 16, for which 
saturation effects do not play a significant role. 

The first thing that we notice is that the dynami- 
cal and valence contributions are roughly additive, in 
the sense that their sum gives back to full signal, in 
the range of fields shown in the figure. The additivity, 
which is expected in the limit of small fields (see dis- 
cussion in Sec. Ill CI) , is verified within errors for b < 8 
(\e\B < (500 MeV) 2 ), while small deviations appear be- 
yond. Notice that this threshold coincides with that 
above which a purely quadratic fit for r(B) does not 
work (see next subsection) and quartic terms in B be- 
come important, in agreement with the argument given 
in Sec. IITCI 

Once clarified that it is sensible, in the explored range 
of fields, to divide magnetic catalysis into a valence and a 
dynamical contribution, from Fig. 2] we learn that the dy- 
namical one is roughly 40% of the total signal, at least for 
the discretization and quark mass spectrum adopted in 
our investigation. That means that numerical studies in 
which the magnetic field is not included in the sampling 
distribution (quenched or partially quenched) may miss a 
substantial part of magnetic catalysis, in a measure larger 
than other systematic effects due to quenching, which are 
typically of the order of 20%. 
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FIG. 4: Relative increment of the average of the u and d quark 
condensates as a function of the magnetic field. We report 
separately r(B), r val {B), r dyn {B) and r val {B) + r dyn (B). 

In Fig. [5] we have also plotted results for the difference 
between the u and the d condensates, which increases as 
a function of B, indicating an increasing breaking of fla- 
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vor symmetry. The fact that the dynamical contribution 
is equal for the two quarks and the approximate additiv- 
ity discussed above implies that such difference should be 
roughly unchanged if we consider just the valence contri- 
bution: that can be verified again from Fig. [5j 
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FIG. 5: Difference of the u and d quark condensates, normal- 
ized by the zero field condensate, as a function of the magnetic 
field and computed respectively on configurations sampled by 
taking or not taking into account the magnetic field in the 
fermion determinant. 



C. Comparison with xPT and model predictions. 

One of the purposes of our investigation is to compare 
our results with various analytic studies based on low 
energy or model approximations of QCD. Most of those 
studies make reference to the average quark condensate 
and not to the u or d condensates separately. Among the 
various existing predictions, one of the first was based 
on the analysis of the Nambu - Jona-Lasino model [l2| 
and predicted a quadratic increase of the condensate as 
a function of the magnetic field, i.e. r(B) oc B 2 . 

A first prediction based on chiral perturbation theory 
has been proposed in Rcf. ff9| 



r{B) 



log(2) \e\B 

\&TT 2 F 2 



(28) 



and is valid only in the chiral limit, i.e. — 0, and for 
\e\B <C Aq CD ; in the limit of strong fields, instead, the 
authors of Ref. [19] have predicted a power law behavior 
r(B) oc \B\ 3 / 2 [19]. Corrections to Eq. ([28]), based on a 
two-loop computations, have been given in Ref. p3| . 

The authors of Ref. [27j have gone beyond the lim- 
itation = 0, presenting a xPT computation which 
is valid for generic values of m 2 ,/{\e\B) 1 even if still for 
\e\B <C Aq CD . The prediction in this case is: 



e LB 



Ih{v) 



log 



i^(log(27r) + ylog(!) 

l + y 



y-21o g r 



(30) 



Notice that Ih{v) — > 1 as y — > 0, i.e. in the chiral limit, 
in agreement with Eq. (l28l) . 

Recently various predictions have been proposed, 
based on the holographic AdS/CFT correspondence (28l . 
l3(il - [3"2l |: the increase in chiral symmetry breaking is con- 
firmed in all cases, with a dependence of the chiral con- 
densate on B which ran ges f rom quadratic (3l| to a power 
law, e.g. r(B) oc \B\ 3 ^ 2 [32j]. 

Existing lattice determinations have reported a linear 
behaviour for SU(2) pure gauge theory [53| and a power 
law behavior r{B) oc B v (with v ~ 1.6) for the SU(3) 
pure gauge theory [H3] • 

Regarding the small field behavior, one can state on 
general grounds that, since by charge conjugation sym- 
metry the chiral condensate must be an even function 
of B, if the theory is analytic at B = then the chiral 
condensate can be written as a Taylor in expansion in 
B 2 , hence for small enough fields the dependence must 
be quadratic. 

That is indeed in agreement with many model predic- 
tions and is not true only in some particular cases: for 
instance the prediction from xPT in Eq. (|28|) 19] is lin- 
ear, since the analyticity requirement is violated in this 
case due to the fact that the pion mass is set to zero. 
Let us consider instead the prediction from Ref. [27j j . 
which is valid for generic values of m 2 /(\e\B). While 
for m 2 /(\e\B) — > Eq. (f29|) gives back a linear behav- 
ior as in Eq. (|28p . in the opposite limit \e\B m 2 , i.e. 
(i << 1), we find, by expanding Eq. (|30)) in powers of 

powers of i, that 



r(B) 



\e\B) 



(31) 



(29) 



in agreement with the general expectation. Such behav- 
ior, quadratic for small fields and linear for larger fields, 
has been found also in a recent study based on the linear 
sigma model [33[ . 

In Fig. [6]we report the relative increment of the u and 
d condensates and of their average in a restricted region 
for which we expect that saturation effects are not impor- 
tant. We remind that in our case the variable represent- 
ing the magnetic field is the dimensionless parameter 6, 
and the conversion to physical units is given by Eq. ([7]) , 
which for our lattice size reads \e\B — (67r/256) b/a 2 , 
with a ~ 0.3 fm, i.e. \e\B ~ b (180 MeV) 2 . 

It is apparent by eye that Eq. badly fits with our 
data, indeed a linear fit gives unacceptable values for 
the x 2 /d.o.f. test (e.g. x 2 / d -°- f - - 190/6 for 1 < b < 7). 
This is expected, since in our case ^ 0. However, data 
at larger values of B show an approximate linear behavior 
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^min 




xVd.o.f 


A B (MeV) 




m in 




X7ndf 


m-n (MeV) 


[MeV) 


l 


4 


0.45 


903(11) 




l 


8 


0.63 


457(59) 


54.7(5.6) 


l 


5 


0.52 


892(6) 




l 


9 


0.93 


392(39) 


61.7(4.2) 


l 


6 


0.67 


899(5) 




l 


10 


0.89 


374(27) 


63.8(2.9) 


l 


7 


0.98 


907(4) 




l 


11 


0.80 


369(20) 


64.3(2.1) 


l 


8 


1.7 


914(4) 




l 


12 


0.79 


359(15) 


65.6(1.6) 


l 


9 


4.6 


923(5) 




l 


13 


0.97 


344(14) 


67.2(1.4) 


l 


10 


8.4 


933(5) 




l 


14 


1.40 


328(14) 


69.1(1.4) 


l 


11 


12 


940(5) 




l 


16 


1.92 


310(13) 


71.1(1.3) 


l 


12 


21 


950(6) 




l 


24 


17.4 


227(22) 


80.8(2.3) 


l 


13 


32 


960(6) 




4 


14 


1.00 


320(12) 


69.9(1.2) 


l 


14 


48 


969(7) 




6 


14 


0.97 


313(13) 


70.6(1.3) 


l 


16 


86 


983(8) 




8 


14 


0.83 


298(15) 


72.2(1.5) 



TABLE II: Results from a quadratic fit r(B) = (\e\B/A%) 2 
to our data. b m i n and 6 max indicate respectively the minimum 
ad maximum values of the magnetic fields included in the fit. 

and if we discard the smallest values of b a function r(b) = 
a + cii b reasonably fits our data: for instance, a fit 
including fe m j n = 6 and 6 max = 16 gives ciq = —0.063(2), 
ai = 0.0195(2) and x 2 /d.o.f. ~ 0.68. 

Next we have tried to check if a quadratic behavior 
r(B) = (|e|B/A|) 2 fits better, at least for small enough 
fields. Results are reported in Table |TTJ For b < 8, i.e. 
\e\B < (500 MeV) 2 , the quadratic fit looks good and 
stable, with A# ~ 900 MeV. That is in agreement with 
the general expectation for the case of small fields dis- 
cussed above (even if 500 MeV is not small with respect 
to Aqcd)- We notice that, in agreement with the ar- 
gument given in Sec. Ill CI the range of validity of the 
quadratic fields roughly coincides with the range in which 
the dynamical and the valence contribution are additive 
(see Sec. HUB]). 
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FIG. 6: Relative increment of the quark condensate as a func- 
tion of the magnetic field. We report separately data for the 
u and d quarks as well as for the average of the two, together 
with our best fit reported in the 6th line of Table IIHI 

We have then tried to fit our data with the predic- 
tion of Ref. [27j, as reported in Eq. (|2"9")l . We have ob- 



TABLE III: Results from a fit of our data to Eq. (J29J) . The 
notation for 6 m m and 6 ma x is as in Table ITT1 



&min 


frmax 


X7ndf 


As (MeV) 


V 


l 


4 


0.23 


784(54) 


2.33(19) 


l 


5 


0.19 


817(25) 


2.23(9) 


l 


6 


0.84 


900(42) 


2.00(12) 


l 


7 


0.92 


940(31) 


1.90(8) 


l 


8 


0.97 


965(24) 


1.85(6) 


l 


9 


1.7 


1006(26) 


1.75(6) 


l 


10 


1.9 


1031(21) 


1.70(5) 


l 


11 


1.9 


1045(18) 


1.67(4) 


l 


12 


2.3 


1064(15) 


1.63(4) 


l 


13 


3.1 


1083(15) 


1.59(4) 


l 


14 


4.3 


1104(16) 


1.55(4) 


l 


16 


6.4 


1131(16) 


1.50(4) 


l 


24 


35.8 


1254(30) 


1.29(5) 



TABLE IV: Results from a fit r(B) = (\e\B/A%) u . The no- 
tation for 6 m in and & ma x is as in Table [TT| 



tained reasonable fits only if both and F n are treated 
as independent free parameters, results are reported in 
Table IIIIl Data are well described by the prediction in 
Eq. (f2"9")) over a wide range of values of \e\B, including 
b = 14 (i.e. \e\B - 700 MeV), even if the fitted values of 
F„ and m w are not very stable as the range is modified. 

Typical values of the fit parameters are ~ 300 — 
400 MeV and F w - 60 - 70 MeV. The fitted pion mass is 
somewhat larger than the value obtained, with the same 
discretization settings, by measuring meson correlators, 
i.e. ~ 200 MeV [42||, however that is quite reasonable 
if we take into account the many discretization system- 
atic effects which affect our simulations. The most rele- 
vant comes from the explicit flavour symmetry breaking 
induced by the staggered discretization: the three pions 
are not degenerate in mass and what is determined by 
meson correlator measurements is just the lowest pion 
mass: it is therefore likely that the xPT prediction still 
holds, but with an heavier effective pion mass. 

Regarding i*V, the fitted values are about ~ 20-30% 
lower than the expected physical value, F„ ~ 93 MeV. 
We notice that F^ enters Eq. (j2"9"]l only in the prefac- 
tor, hence its value is surely affected by the systematic 
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uncertainty in the overall normalization factor for r(B), 
which is of the order of 10%. There are also corrections 
expected from the fact that we are not at zero tempera- 
ture: the authors of Ref. [H| predict F%(T) ~ F 2 -T 2 /6, 
however in our case T ~ 40 MeV and that can account 
for at most a 2% discrepancy from the physical value. 
The non-physical large value of m n can also affect F T , 
but xPT would predict an increased value of F v (65j | . 

There are however many other possible sources of sys- 
tematic uncertainties, including the fact that the xPT 
prediction of Ref. [27| has been obtained in the low en- 
ergy limit \e\B -C Aq CD , a condition which is violated 
in our explored range of fields. We have verified that 
other two-parameter functions, which allow to fix inde- 
pendently the curvature at B = and the asymptotic lin- 
ear behavior for larger fields (like Eq. (|29p when F n and 
ttt-tt are treated as independent parameters) work equally 
well. For instance the function 



r(b) = cob atan(ci&) 



(32) 



fits well in the whole range of explored fields, 6 m j n = 1 
and b max = 16, with c = 0.0136(2), ei = 0.140(5) and 
X 2 /d.o.f. ~ 0.77. 

To compare with the analysis performed in Ref. [5~j| . 
we have also investigated if a power law behaviour r(B) — 
(\e\B/K 2 B ) v can fit our data. Results are reported in 
Table ITVl Reasonable fits are obtained only for a range 
of fields including 6 = 8. That coincides more or less with 
the range for which also the quadratic fit works well, and 
indeed values obtained for v in this range are roughly 
compatible with v = 2. 

We are not able to say much about the strong field 
regime, \e\B ^> Aq CD , since saturation effects make such 
regime unaccessible to our present investigation: much 
smaller lattice spacings should be used to that aim. 



D. Density of zero modes 

Let us finally discuss our results for the densities of 
zero modes of the Dirac operator, p u (0) and p<j(0), ob- 
tained for the u and d quarks respectively. We have de- 
termined such densities following Eq. (|20[) . i.e. deter- 
mining on our ensemble of configurations, sampled with 
a dynamical quark mass am = 0.01335, the average of 
Ti(M~ 1 (am')) for various different values of m! and then 
extrapolating to am' — 0. In particular we have chosen 
am' = 0.007, 0.01335, 0.02 and 0.03. A quadratic extrap- 
olation in ml works well in all cases. 

In Table [V] we report data for the relative increments 
f u and fd of jO„(0) and Pd(fi) as a function of B, also 
for the cases in which only the valence or the dynamical 
contributions are taken into account; the average quanti- 
ties are plotted in Fig. [7] We notice, comparing Table IVl 
with TablcU that the relative increment of the zero mode 
density, f, is generally larger than the corresponding in- 
crement in the condensate, r. Moreover Fig. [7] shows 
that, similarly to what happens for the condensate, it is 



b 


r„(6) 


fdQ>) 


fT\b) 


f v d al {b) 




1 


0.003(4) 


0.000(4) 


0.002(3) 


0.002(3) 


0.002(3) 


2 


0.013(4) 


0.002(5) 


0.007(3) 


0.003(3) 


0.003(3) 


3 


0.026(4) 


0.014(5) 


0.022(4) 


0.007(4) 


0.009(4) 


4 


0.042(3) 


0.024(4) 


0.034(4) 


0.010(5) 


0.012(4) 


5 


0.069(3) 


0.034(5) 


0.047(3) 


0.016(3) 


0.018(3) 


6 


0.080(3) 


0.040(7) 


0.068(3) 


0.021(4) 


0.027(4) 


7 


0.117(3) 


0.062(4) 


0.087(3) 


0.031(4) 


0.036(3) 


8 


0.144(5) 


0.076(5) 


0.113(6) 


0.037(5) 


0.045(5) 


9 


0.173(5) 


0.088(4) 


0.138(5) 


0.044(4) 


0.050(3) 


10 


0.201(3) 


0.101(7) 


0.158(5) 


0.052(5) 


0.056(3) 


11 


0.236(5) 


0.124(4) 


0.185(5) 


0.063(5) 


0.069(4) 


12 


0.264(4) 


0.142(5) 


0.207(4) 


0.074(5) 


0.076(3) 


14 


0.319(5) 


0.178(8) 


0.256(4) 


0.091(4) 


0.097(3) 


16 


0.372(5) 


0.211(5) 


0.302(6) 


0.115(4) 


0.110(3) 



TABLE V: Relative increment of the density of zero modes 
of the Dirac operator for the u and d quarks and for various 
magnetic field values. We report full data, as well as valence 
and dynamical contributions separately. 



possible to make a separation of the total increment into 
a valence and a dynamical contribution, which sum ap- 
proximately to the total increment. Also in this case the 
change in the distribution of the non-Abelian gauge fields 
(dynamical contribution) accounts for about 30-40% of 
the total increase. 
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FIG. 7: Relative increment of the average of the p u (0) and 
Pd(0) as a function of the magnetic field. We report separately 
r{B),f val (B),f dyn (B) and f val {B) +? dyn (B). 



IV. CONCLUSIONS 

In this study we have approached the issue of mag- 
netic catalysis by numerical lattice simulation of Nj = 2 
QCD. We have adopted a rooted standard staggered dis- 
cretization of the fermion action and a plaquette pure 
gauge action on a symmetric 16 4 lattice, corresponding 
to a lattice spacing a ~ 0.3 fm, a (Goldstone) pion mass 
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m v ~ 200 MeV and a temperature well below the decon- 
finement/chiral restoring one (T ~ 40 MeV). We have 
studied the breaking of chiral symmetry as a function of 
a constant and uniform magnetic field directed along the 
z axis. Explored magnetic fields, allowed by the toroidal 
geometry and for which distortion (saturation) effects 
due to the lattice discretization are not significant, range 
from \e\B ~ (180 MeV) 2 to \e\B ~ (700 MeV) 2 . 

We have shown that, in the range of explored fields, 
it is possible to divide magnetic catalysis into a con- 
tribution coming from the modified distribution of non- 
Abelian gauge fields, induced by dynamical quark loop 
effects, that we have called dynamical contribution, and a 
valence contribution, determined by measuring the con- 
densate on gauge configurations sampled with the un- 
modified distribution. The first term, which is missed 
by quenched or partially quenched studies, accounts for 
about 40% of the total increase in the quark condensate. 
Results obtained for the density of zero modes looks quite 
similar. 

Regarding the dependence of the condensate on the 
magnetic field, we have shown that a quadratic behavior, 
which is expected in the limit of small magnetic fields, 
describes well our data for \ e\B up to ~ (500 MeV) 2 . 
The xPT prediction of Ref. [27[ fits data over a wider 
range, but only if the pion decay constant is treated as 
an independent free parameter. 



Our investigation can be improved in several respects. 
Lattice artifacts may affect our results in various differ- 
ent ways, ranging from the presence of a distorted me- 
son spectrum in the adopted rooted staggered fermion 
formulation, to residual renormalization effects and pos- 
sible residual saturation effects in the explored range of 
fields. An improved lattice formulation and a finer spac- 
ing a would allow to check for such artifacts, to test the 
correct scaling to the continuum limit of our results and 
to explore larger values of the magnetic field. A larger 
spatial volume would instead allow for a finer quantiza- 
tion of \e\B and a better investigation of the small field 
region. It would be also interesting to explore different 
choices of the quark mass spectrum, in order to see how 
the separation of magnetic catalysis into a dynamical and 
a valence contribution depends on the dynamical quark 
masses. We plan to address such issues in future studies. 
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